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FIG. 1. Phase diagram of the CG in the u — T plane. 

ary terms, it is easy to computciS the inverse dielectric 
tensor , which is the equivalent of the helicity modu- 
lus of the corresponding XY model 00 

( -.=2»(v»(^)) a ,„ 

^vr^w^Y), ( 2) 



T \ \ L J \ L 

where V and V" are first and second derivatives of V. 
For an isotropic system, t7l = (r l 5^ v . According to 
the KT instability criterion,B in the insulating phase we 
must have the inequality e _1 > 4T. This will enable us to 
set an upper bound on the transition to the conducting 
phase. 

We carry out standard Metropolis Monte Carlo sim- 
ulations for various values of u and T. At each point 
an initial 20,000 MC passes are discarded to equilibrate 
with an additional 128,000 MC passes to compute aver- 
ages. Errors are estimated from block averages. 

The u —JT phase diagram, as found previously by Lee 
and TeitelJ3 is shown in Fig.[lj. At low u, upon increasing 
T, there is a KT transition to a conducting liquid. At 
low T, upon increasing u, there is a first order transition 
to an insulating charge solid at uq — tt/8. Increasing 
T within this charge solid gives first a KT transition to 
a conducting solid, followed by an Ising melting of the 
solid. The two KT lines meet the first order line at the 
same T* ~ 0.125, slightly lower than the tricritical point 
where the first order and Ising lines meet. 

In Fig. we present our new results for this phase 
diagram in the p — T plane. The coexistence boundary is 
determined as follows. The low temperature branches are 
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FIG. 4. Inverse dielectric function e _1 (T) for various lattice 
sizes L, for u = 0.39 just below uo = tt/8. The intersection 
with the line AT gives the KT upper bound on the transition 
temperature T* . 



was observed in e _1 . 

As a more precise test that the transition is indeed 
KT-like, we use the-Jjuaite size scaling procedure of We- 
ber and MinnhagenllSEJ. Precisely at the KT transition 
temperature, the finite size dependence of e _1 is given 
by, 



- 1 (L) = e-\oc) 



1 

1 



21nL 



(3) 



Fitting data for various L at fixed T to Eq.([|), with 
e _1 (oo) and c as free parameters, one determines Tkt as 
that temperature which gives the smallest x 2 error of the 
fit. The fitted value of e _1 (oo) at the Tkt so determined, 
should then turn out to be precisely e _1 (oo) = 4Tkt, so 
as to obey the universal KT prediction. We show the 
results of such a fit below. In Fig.||a we show the x 2 of 
the fit versus temperature, where we have included sizes 
L = 8 — 64, 12 — 64, and 16 — 64 in the fit. In each case, 
the minimum of \ 2 occurs at T* — 0.1235. In Fig.|j^6 we 
show the corresponding fitted value of e _1 (oo)/T versus 
temperature. We see that e _1 (oo)/T = 4 at essentially 
the same T* where \ 2 has its minimum. This analysis 
strongly supports the transition as being of the KT type. 

Finally, in Fig. ^ we show histograms of the density p 
for u = 0.39 and several values of T passing through T* , 
for the largest system size we have studied, L = 64. We 
have chosen the values of T shown in Fig.|| so that the 
histograms between neighboring values have a significant 
overlap. As is clearly seen, all the distributions are single 
peaked. There is no sign at all of the bimodal distribution 
that would characterize a first order transition. While we 
cannot rule out the possibility of a weak first order tran- 
sition with a £(T*) > 64, our results clearly suggest that 
the "KT" line at small u does indeed remain a second or- 
der Kosterlitz-Thouless transition all the way up to uq. 



Our lattice simulation has the advantage that it is eas- 
ier to equilibrate at high densities than continuum sim- 
ulations. However the presence of the discrete lattice 
does have a strong influence on the location of the phase 
boundaries. In particular, the coexistence boundary of 
Fig. H at small T < 0.1 is well described by a simple model 



FIG. 6. Histograms of charge density p for u = 0.39 just 
below uo = 7r/8. Histograms remain single peaked as one 
passes though the transition at T* ~ 0.1235 



of excited isolated dipoles and quadrupoles. For the 
lower branch we get p-{T) = ^(2N d e-^ Ed +4N q e~l 3E i), 
where Ed and E q are the excitation energies of an iso- 
lated dipole and quadrupole respectively, and Nd = 4L 2 
and N q = 21? are the number of ways of putting them 
down on the square lattice. For the upper branch we get 
p+(T) = 1 - ^{2N d e~ 0E '" +AN' q e- 0E 'i) where E' d and 
E' are the energies for removing an isolated dipole and 
quadrupole, and N' d = 21? and N' q = L 2 are the number 
of ways this may be done. Clearly these expressions will 
change with the geometry of the discretizing lattice, or if 
a continuum is used. Nevertheless, at the very low den- 
sities p ~ 0.004 and high temperatures T ~ 0.25 where 
Levin et al. estimate a tricritical point, we would be very 
surprised if the lattice is qualitatively different from the 
continuum. 

Using a discrete lattice also has the effect that it tends 
to stabilize the charge solid phase above uq to high tem- 
peratures. Indeed, our charge solid only melts after it has 
already become conducting via a KT transition arising 
from the excitation and diffusion of vacancies through- 
out the solid. The present model does not possess any 
first order transition from an insulating tn a conduct- 
ing phase. In contrast, recent simulationsEallj of the CG 
in a flat continuum with periodic boundary conditions 
find that the charge solid phase is melted at any finite 
temperature. Here, a first order line separates the in- 
sulating gas from a dense conducting liquid, ending at 
a critical point at relatively low temperature and high 
density: (T c , p c ) = (0.056, 0.21) according to Ref. ||, 
and T c — 0.032 according to Ref. [l9]. Similar results 
were found earlier for the continuum CG on the surface 
of a sphered (T c ,p c ) = (0.087,0.11). In these models, 
the KT line ends either at, or near this critical point. 

Although the geometry of the the CG system clearly 
affects the location of the end of the KT transition line, 
it is interesting to compare our results for the square lat- 
tice with the predictions of the continuum self-consistent 
screening theory of Mipjahagen and Wallin.El To do so, 
it is necessary to noteoD that if the interaction G{f) 
on the lattice is chosen so as to asymptotically match 
the continuum — lnr as r — > oo, then the lattice CG 
with u = acts like a continuum model with a chemical 
potential p = -Aln(8e 2 T) ~ 0.8085. Thus a chemi- 
cal potential uq — n/8 on the grid acts like a chemi- 
cal potential p = uq — po = —0.416 in the continuum. 
Minnhagen and Wallin predict that the KT line will end 
at T* = 0.144 and fugacity z* = 0.054, giving a chemical 
potential p = T* \nz* = -0.420. 

To conclude, we show that the 2D neutral CG contin- 
ues to have an ordinary KT transition up to high densi- 
ties p ~ 0.1, in contrast with recent theoretical estimates. 
This result, obtained here for the square lattice CG, is 
consistent with recent simulations in the continuum. 
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